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Abstract 

We reexamine a very classical problem, the spinning behavior of the 
tippe top on a horizontal table. The analysis is made for an eccentric 
sphere version of the tippe top, assuming a modified Coulomb law for 
the sliding friction, which is a continuous function of the slip velocity Vp 
at the point of contact and vanishes at Vp = 0. We study the relevance of 
the gyroscopic balance condition (GBC), which was discovered to hold 
for a rapidly spinning hard-boiled egg by Moffatt and Shimomura, to the 
inversion phenomenon of the tippe top. We introduce a variable £ so that 
£ = corresponds to the GBC and analyze the behavior of £. Contrary 
to the case of the spinning egg, the GBC for the tippe top is not fulfilled 
initially. But we find from simulation that for those tippe tops which will 
turn over, the GBC will soon be satisfied approximately. It is shown that 
the GBC and the geometry lead to the classification of tippe tops into 
three groups: The tippe tops of Group I never flip over however large a 
spin they are given. Those of Group II show a complete inversion and 
the tippe tops of Group III tend to turn over up to a certain inclination 
angle Of such that 9f<n, when they are spun sufficiently rapidly. There 
exist three steady states for the spinning motion of the tippe top. Giving 
a new criterion for stability, we examine the stability of these states in 
terms of the initial spin velocity no- And we obtain a critical value n c 
of the initial spin which is required for the tippe top of Group II to flip 
over up to the completely inverted position. 
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1 Introduction 



Spinning objects have historically been interesting subjects to study. The spin 
reversal of the rattleback pQ (also called a celt or wobblestone) and the behavior of 
the tippe top are typical examples. In the latter case, when a truncated sphere with 
a cylindrical stem, a so-called 'tippe top', is spun sufficiently rapidly on a table with 
its stem up, it will flip over and rotate on its stem. This inversion phenomenon has 
fascinated physicists and has been studied for over a century OlSllllElEllZllHlini- 

In the present paper we revisit and study this very classical problem from a 
different perspective. Recently the riddle of spinning eggs has been resolved by 
Moffatt and Shimomura [MS] [10J. They discovered that if an axisymmetric body, 
such as a hard-boiled egg, is spun sufficiently rapidly, a 1 gyroscopic balance' condition 
(GBC) holds and that under this condition the governing equations of the system 
are much simplified. In particular, they derived a first-order ordinary differential 
equation (ODE) for 9, the angle between the axis of symmetry and the vertical axis, 
and showed for the case of a prolate spheroid that the axis of symmetry indeed 
rises from the horizontal to the vertical. Then the spinning behavior of egg-shaped 
axisymmetric bodies, whose cross sections are described by several models of oval 
curves, was studied under the GBC by one of the present authors jTTJ. 

The tippe top is also an axisymmetric body and shows the similar behavior as 
the spinning egg. Then one may ask: does the GBC also hold for the tippe top? 
If so, how is it related to the inversion phenomenon of the tippe top? In the first 
half of this paper we analyze the spinning motion of the tippe top in terms of the 
GBC. Actually the GBC is not satisfied initially for the tippe top, contrary to the 
case of the spinning egg. The difference comes from how we start to spin the object: 
we spin the tippe top with its stem up, in other words, with its symmetry axis 
vertical while the egg is spun with its symmetry axis horizontal. In this paper we 
perform our analysis taking an eccentric sphere version of the tippe top instead of a 
commercially available one, a truncated sphere with a cylindrical stem. In order to 
examine the GBC of the tippe top more closely, we introduce a variable £ so that 
£ = corresponds to the GBC, and study the behavior of £. Numerical analysis 
shows that for the tippe tops which will turn over, the variable £, starting from a 
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large positive value £o> soon takes negative values and fluctuates around a negative 
but small value £ m such that |£ m /£o| ~ 0. Thus for these tippe tops, the GBC, which 
is not satisfied initially, will soon be realized but approximately. On the other hand, 
in the case of the tippe tops which will not turn over, £ remains positive around 
£o or changes from positive £o to negative values and then back to positive values 
close to £o again. We find that the behavior of £ is closely related to the inversion 
phenomenon of the tippe top. Once £ fluctuates around the value £ m , the system 
becomes unstable and starts to turn over. 

Under the GBC the governing equations for the tippe top are much simplified 
and we obtain a first-order ODE for 9, which has the same form as the one derived 
by MS for the spinning egg. Then, this equation for 9 and the geometry lead to the 
classification of tippe tops into three groups, depending on the values of ^ and j|, 
where A and C are two principal moments of inertia, and a is the distance from the 
center of sphere to the center of mass and R is the radius of sphere. The tippe tops 
of Group I never flip over however large a spin they are given. Those of Group II 
show a complete inversion and the tippe tops of Group III tend to turn over up to a 
certain inclination angle 9f such that 6j <tt, when they are spun sufficiently rapidly. 
This classification of tippe tops into three groups and its classificatory criteria totally 
coincide with those obtained by Hugenholtz P] and Leutwyler p], both of whom 
resorted to completely different arguments and methods. 

In the latter half of this paper we study the steady states for spinning motion of 
the tippe top and examine their stability (or instability). It is well understood that 
the main source for the tippe top inversion is sliding friction jSEj, which depends 
on the slip velocity Vp of the contact point between the tippe top and a table. Often 
used is Coulomb friction (see Eq. (j2.1Hj) ). In fact, Coulomb friction is practical when 
\vp\ is away from zero, but it is undefined for vp = 0. However, we learn that at 
the steady state of the tippe top, the slip velocity v P necessarily vanishes. In order 
to facilitate a linear stability analysis of steady states and also to study the motion 
of the tippe top as realistically as possible, we adopt in our analysis a modified 
version of Coulomb friction (see Eq. ()2.14)l ). which is continuous in up and vanishes 
at v P = 0. 
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Actually the steady states of the tippe top and their stability were analyzed by 
Ebenfeld and Scheck [ES] [7], who assumed a similar frictional force which is con- 
tinuous at vp = 0. They used the total energy of the spinning top as a Liapunov 
function. The steady states were found as solutions of constant energy. And the 
stability or instability of these states were judged by examining whether the Lia- 
punov function assumes a minimum or a maximum at these states. Also recently, 
Bou-Rabee, Marsden and Romero [BMR] [§| analyzed the tippe top inversion as 
a dissipation-induced instability and, using the modified Maxwell-Bloch equations 
and an energy- momentum argument, they gave criteria for the stability of the non- 
inverted and inverted states of the tippe top. 

We take a different approach to this problem. First, in order to find the steady 
states for spinning motion of the tippe top, we follow the method used by Moffatt, 
Shimomura and Branicki [MSB] for the case of spinning spheroids [1.2] • Then the 
stability of these steady states is examined as follows: Once a steady state is known, 
the system is perturbed around the steady state. Particularly we focus our attention 
on the variable 8, which is perturbed to 8 = 8 s + 58, where 8 S is a value at the steady 
state and 58 is a small quantity. Using the equations of motion, we obtain, under the 
linear approximation, a first-order ODE for 58 of the form, 58 = H s 58, where H s is 
expressed by the values of dynamical variables at the steady state. Thus the change 
of 58 is governed by the sign of H s . If H s is positive (negative), \58\ will increase 
(decrease) with time. Therefore, we conclude that when H s is negative (positive), 
then the state is stable (unstable). Using this new and rather intuitive criterion we 
argue about the stability of the steady states in terms of the initial spin velocity tiq 
given at the position near # = 0. We observe that our results on the stability of the 
steady states are consistent with ones obtained by ES and MSB. Then we obtain a 
critical value n c of the initial spin which is required for the tippe top of Group II 
to flip over up to the completely inverted position at 8 = n. Finally we confirm by 
simulation our results on the relation between the initial spin no and the stability 
of the steady states. 

The paper is organized as follows: In Sec. 2 we explain the notation and geometry 
used in this paper, and give all the necessary equations for the analysis of the 
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spinning motion of the tippe top. In Sec. 3 we discuss about the GBC and its 
relevance to the inversion phenomenon of the tippe top. We also show that the 
assumption of the GBC leads to the classification of tippe tops into three groups. 
Then in Sec. 4 we study the steady states for the spinning motion of the tippe 
top and examine their stability. Sec. 5 is devoted to a summary and discussion. 
In addition, we present four appendices. In Appendix A, the equations of motion 
which are used to analyze the spinning motion of the tippe top are enumerated. In 
Appendix B, it is shown that intermediate steady states for the tippe tops of Group 
II and Group III are stable when an initial spin n{9 = 0) falls in a certain range. 
In Appendix C we demonstrate that our stability criterion for the steady state is 
equivalent to the one obtained by ES. And finally, in Appendix D, we show that 
our results on the stability of the vertical spin states are consistent with the criteria 
derived by BMR. 



2 Equations of motion for tippe tops 

A commercially available tippe top is usually a truncated sphere with a cylindrical 
stem. Instead we perform our analysis taking a loaded (eccentric) sphere version 
of the tippe top. The center of mass is off center by a distance a. There are no 
qualitative differences between the two. But if applied to the case of a commercial 
tippe top with a stem, our assertions would be valid up to the point when the stem 
touched the table surface. 

Fig. ^ shows the geometry. An axisymmetric tippe top spins on a horizontal 
table with point of contact P. We will work in a rotating frame of reference OXYZ, 
where the center of mass is at the origin, O. The center S of the sphere with radius 
R is at a distance a from the origin. The symmetry axis of the tippe top, Oz, and the 
vertical axis, OZ, define a plane II, which precesses about OZ with angular velocity 
Q(t) = (0, 0, Q). Let (0, 9, ip) be the Euler angles of the body relative to OZ. Then 
we have Q = <fi, where the dot represents differentiation with respect to time, and 9 
is the angle between OZ and Qz. We choose the horizontal axis OX in the plane II 
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Figure 1: A loaded sphere (eccentric) version of the tippe top. The center of mass 
O is off center (S) by distance a. The tippe top spins on a horizontal table with 
point of contact P. Its axis of symmetry, Qz, and the vertical axis, OZ, define a 
plane II, which precesses about OZ with angular velocity f2(t) = (0, 0, ft). OXYZ is 
a rotating frame of reference with OX horizontal in the plane II. The height of O 
above the table is h{6) = R — acos9, where R is the radius. The position vector of 
P from O is X P = (X P , 0, Z P ), where X P = § and Z P = -h{6). 

and thus OY is vertical to II and inward. 

In a rotating frame of reference Oxyz, where Ox is in the plane II and perpendic- 
ular to the symmetry axis Qz and where Oy coincides with OY, the tippe top spins 
about Qz with the rate ip. Since f2 is expressed as f2 = — VLsmOx + VtcosOz 
in the frame Oxyz, the angular velocity of the tippe top, uj, is given by uj = 
— Qsin8x + 9y + nz. Here x, y, and z are unit vectors along Ox, Oy, and Oz, 
respectively, n(t) is given by n = Qcosd + ip. The Ox and Oy are not body-fixed 
axes but are principal axes, so that the angular momentum, L, is expressed by 
L = —AQ sin 6x + A6y + Cnz, where (A, A, C) are the principal moments of inertia 
at O. Using the perpendicular axis theorem and the parallel axis theorem, we see 
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that A/C > | for any axisymmetric density distribution. 

The coordinate system Oxyz is obtained from the frame OXYZ by rotating the 
latter about the OY (Oy) axis through the angle 9. Hence, in the rotating frame 
OXYZ, u) and L have components 



(^(n-ncos9)sin9,9 , sin 2 # + n cos , (2.1) 
L = ((Cn- AVLcosO) sinO, AO , Ansin 2 9 + Cncos9^j , (2.2) 

respectively. The evolution of L is governed by Euler's equation 

dL 

— + QxL = X P x(N + F), (2.3) 

where Xp is the position vector of the contact point P from O, N is the normal 
reaction at P, N = (0,0, N), with iV being of order Mg, the weight, and F = 
(Fx,Fy,0) is the frictional force at P. We consider only the situation in which the 
tippe top is always in contact with the table throughout the motion. Since the point 
P lies in the plane II, Xp has components (Xp, 0, Zp), which are given by 

Z P = -(R- a cos 6) = -h(6) , (2.4a) 
X P = a sm9 =^-, (2.4b) 

de 

where h(6) is the height of O above the table. The components of ()2.3j) are expressed, 
respectively, as 

L x - ttL Y = h(e)F Y , (2.5a) 
Ly + VLLx = -asineN -h(6)F x , (2.5b) 
L z = asindFy . (2.5c) 

In terms of 9, Q, and n the above equations are rewritten as 

AQ sin 9 = (Cn - 2 AQ cos 9)9 + (a- R cos 9) F Y , (2.6a) 
A9 = -n(Cn- Ancos9)sin9 -a sin 9N -h(9)F x , (2.6b) 
Ch = Rsin9F Y . (2.6c) 



6 



Now it is easily seen from (|2.2|) . (J2.5a|) and (|2.5cj) that there exists an exact 
constant of motion, 

J = -L ■ X P = -L x ^- + L z h(9) (a constant), (2.7) 

d9 

which is valid irrespective of the reaction force (N + F) at the contact point P, in 
other words, whether or not slipping occurs. This so-called "Jellett's constant" [13J 
is typical for the tippe top whose portion of the surface in contact with the table is 
spherical. 

The velocity, v TOt p, of the contact point P with respect to the center of mass O 
is given by v TOtP = u> x X P , and thus has components, 

Vrotpx = -h{6)0 , (2.8a) 
VrotPY = {R(n — Q cos 9) + aQ} sin 9 , (2.8b) 
Vrotpz = —a sin 9 9 . (2.8c) 

The center of mass O is not stationary. Let uq = {uox,uoy,uoz) represent the 
velocity of O, then the slip velocity of the contact point P, v P = (vpx,vpy,vpz), 
is 

Vp = U + ^rotP • (2.9) 

Since uoz = % = —v m tpz , we have v P z = as was expected. 
The equation of motion for the center of mass O is given by 

M^^ + nxuo^j =N + F + W , (2.10) 

where M is the mass of the tippe top and W = (0, 0, —Mg) is the force of gravity. 
In components, Eq. ()2.10j) reads 

M{u ox -nu OY ) = F x , (2.11a) 

M(u Y + ttuox) = F Y , (2.11b) 

Miioz = N-Mg. (2.11c) 

Since ii z = §S F;q. (j2.11cj) gives 

N = M^g + a (9 2 cos9 + 9 sin^j , (2.12) 



which shows that the normal force N is of order Mg when a9 2 , a\9\ <C g. 

We need an information on the frictional force F. It is well understood that the 
sliding friction is the main source for the tippe top inversion [3 E] • So we will ignore 
other possible frictions, such as, rolling friction ^3] and rotational friction which is 
due to pure rotation about a vertical axis . 

Concerning the sliding friction, often used is a Coulomb law, which states that 

F c = -f^N^j . (2.13) 
\ V P\ 

where \i is a coefficient of friction. Another possibility is a viscous friction law, 
which states that the friction is linearly related to vp. Coulomb friction is practical 
when \vp\ is away from zero but it is undefined at Vp = 0. The slip velocity of the 
contact point P necessarily vanishes at the steady state of the tippe top. In order 
to study the motion of the tippe top as realistically as possible and also to facilitate 
a linear stability analysis of steady states, we modify the expression of Coulomb 
friction (|2.1Hj) as 

F = -fiN-^— with MA)| = ^x+^ + A 2 , (2-14) 

so that F is continuous in and vanishes at vp = 0. Here we choose A as a 
sufficiently small number with dimensions of velocity. Note that vpz = and thus 
the Z-component of F is 0. 

This completes the presentation of all the necessary equations for the analysis of 
the motion of tippe tops. We enumerate all these equations in Appendix A. We need 
further the initial conditions. When we play with a tippe top, we usually give it a 
rapid spin with its axis of symmetry nearly vertical. So let us choose the following 
initial conditions for 9 and other angular velocities: 

9 Q = 6(t = Q) small, 9 = 6(t = 0) = 
tt = tt(t = 0) = 0, (2.15) 
tpo = ip(t = 0) large . 

We take 9 = 0.01 ~ 0.1 rad and V'o = 10~150 rad/sec. Recall that the spin n(t) is 
given by n = Q cos 9 + ip, and thus we have hq = n(t = 0) = 10 ~ 150 rad/sec. As for 
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the initial condition for the velocity of the center of mass O, we take 

u = u o {t = 0) = , (2.16) 

since we usually do not give a large translational motion to the tippe top at the 
beginning. 

With the above initial conditions (j2.15|) and (j2.16|) . we analyze the behaviors of 
the tippe top using three angular (|2.6aH2.6cj) and three translational (j2.11aH2.11c|) 
equations of motion, together with the knowledge of the frictional force, a modified 
version of the Coulomb law (j2.14j) . and the velocities (|2.8aH2.8c|) and ()2.9|) . When 
we perform simulations we use the adaptive Runge-Kutta method. 



3 Gyroscopic balance condition 

3.1 The variable f 

We define a variable £ as 

£ = Cn- Attcos6 . (3.1) 

In terms of £, the X- and Z- components of L in ()2.2j) and Jellett's constant J, 
(J2.7j) . are expressed, respectively, as 

L x = £ sin#, L z = f cosfl + Att , (3.2) 
J = -£a sin 2 6 + L z h{6) . (3.3) 

The condition £ = has been introduced by MS ^U] in their analysis of spinning 
hard-boiled eggs, and referred to as the GBC They discovered that the GBC, £ = 0, 
is approximately satisfied for the spinning egg and, using this GBC, they resolved a 
long standing riddle: when a hard-boiled egg is spun sufficiently rapidly on a table 
with its axis of symmetry horizontal, the axis will rise from the horizontal to the 
vertical. We outline how MS found the GBC for the spinning egg ^Hj- The system 
of the spinning egg obeys essentially the same equations of motion as the case of 
the tippe top, to be specific, Eqs. ()2.3|) and (|2.10j) . The F-component of ()2.3|) 
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for the spinning egg is given by (j2.6b|) . with the factor, asin^, being replaced by 
Xp. Because the secular change of 9 is slow and thus \9\ <C Q 2 , the term AO can 
be neglected. Furthermore, in a situation where Q 2 is sufficiently large so that the 
terms involving Q in (j2.6b|) dominate the terms —XpN and —h(6)Fx, Eq. (|2.6b|) is 
reduced, in leading order, to (Cn — AQ cos 9)Q sin 9 = 0. Hence, for sin# ^ 0, we 
arrive at the condition £ = Cn — AQ cos 6 = 0. 

The tippe top shows the similar behavior as the spinning egg. Then one may ask: 
does the GBC also hold for the tippe top? We will show that the answer is "partly 
no" and "partly yes". "Partly no" means that the GBC is not satisfied initially. 
Tippe tops are usually spun with 9 w 0, Q ~ 0, and large ip and, therefore, 
n pa ipQ is large, from which we find that £o = £(^ = 0) ~ Cn Q is large 1 . Thus the 
GBC does not hold at the beginning. However, we will see later that the GBC does 
approximately hold whenever the tippe top rises, which is the meaning of "partly 
yes". In fact, the argument of MS to derive the GBC for the spinning egg can also 
be applied to the tippe top. Thus in a situation where fl is sufficiently large and for 
sin 9 7^ 0, the GBC is expected to be satisfied. On the other hand, in the case of 
the spinning egg, the GBC is approximately satisfied initially. We start to spin an 
egg with its symmetry axis horizontal, that is, with 6*o ~ f, ~ and large Qq. 
Hence we find no ~ and £o ~ for the spinning egg. 

We emphasize that the variable £ initially takes a large positive value for the 
tippe top. But our numerical analysis will show that when a tippe top turns over, 
£ soon makes a rapid transition from large positive values to negative values and 
starts to oscillate about a small negative value. 

Before proceeding with a discussion of this transition of £, let us consider the 
consequences when the GBC is exactly satisfied for the tippe top. 

3.2 Consequences of the exact GBC 

In a situation where Q is sufficiently large and 9 is not in the vicinity of or tt, the 
GBC is realized for the tippe top. Let us consider the case that the exact GBC, 

this paper we always take the initial spin velocity ipo about Oz to be positive and, therefore, 
£o is positive. 
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£ = 0, is satisfied for the tippe top. Then, we have 

J = L z h(6) , (3.4) 

from ()3.3|) . and L z = AQ from the second equation in ([3.2)1 . If the angular velocity 
Q around the vertical axis is reduced and, therefore, L z decreases, Eq. ([3.4)1 tells 
us that the height h{9) of the center of mass from the table increases since J is a 
constant, which means the turning over of the tippe top. Differentiating both sides 
of ()3.4j) by time and using (j2.4b|) and ([2.5cjl . we obtain a first-order ODE for 9, 

J0 = -F Y h 2 (6) . (3.5) 

We assume also that the ^-component of uo, the translational velocity of the 
center of mass O, in ([2.9)) is negligible in the first approximation as compared with 
that of f r otp, and we set vpy = v TOt py. We see that numerical simulation supports 
this assumption. Then, one can use Eq. (j2.8b|) and the GBC to eliminate n and Q, 
and obtain vpy as only a function of the dynamical variable 9 as follows: 

Jsin^ f [A \ J . . 

Since the frictional force F Y is proportional to vpy, we obtain from Eqs. ([3.5H3.6)) . 

6 oc vpy (3.7) 
with a positive proportional coefficient and 

vp Y = sinflja + i?^ - l) cos #j . (3.8) 

Equation ()3.7)) implies that the change of 6 is governed by the sign of vpy. If vpy 
is positive (negative), then 9 will increase (decrease) with time. Therefore a close 
examination of the behavior of vpy as a function of 9 will be important 2 . 

We observe from ()3.8[) that vpy = at 9 = and ir, since sin (9 = at these angles. 
Moreover, vpy may vanish at an other angle, which is given by solving 

(A \ 

a + R[g - lj cos# = . (3.9) 



2 A resemblance of (|3.7f) to a renormalization group equation which appears in quantum field 
theories for critical phenomena and high energy physics is emphasized in Sec. 5. 
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(a) 



Inclination Angle 



Inclination Angle 6 

(b) 




Figure 2: Vpx as a function of 9 for tippe tops of (a) Group I with j| = 0.1 and 
4 = 0.8; (b) Group II with | = 0.1 and 4 = 1; (c) Group III with | = 0.1 and 



Equation ()3.9|) has a solution for if ^ < 1 — -| or 1 + — < ^ and no solution 
otherwise. Accordingly, tippe tops are classified into three groups, depending on the 
values of ^ and -|: Group I with ^ < 1 — j|; Group II with 1 — j| < ^ < 1+J|; and 
Group III with 1 + | < 4. 

We now examine the behaviors of tippe tops belonging to each group, 
(i) Group I (£<1-|) 

Imagine that a familiar top consisting of a circular disk and a stem is located inside 
of a hollow massless sphere. The stem is along the diameter of the sphere whose 
center does not coincide with the center of mass. This toy may belong to Group I. 
Figure El (a) shows a typical graph of vpy for a tippe top of Group I. The graph 
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crosses the line Vp = at an angle 



7T 

R(l-f,)l B '2 



# c = cos" 1 ( — ] and < C < - , (3.10) 



and vpy is positive for C < < tt but negative for < < C . So, the angle 9 C is a 
critical point. If a tippe top of Group I is spun on a table with ipo m and sufficiently 
large f2 and with the initial angle 0q > @c-, then 9 will increase to 7T, which means 
that the body will eventually spin at 9 = n. In the case 9q < 9 C , we will see that the 
body spins at 9 = 0. Depending on the initial value 9q the body will spin at the end 
point 9 = or n. Both ends are stable points. Usually we give a spin to the tippe 
top at a position with 9 rs 0. Spun at 0q ~ 0, the tippe top of Group I does not 
turn over however large a spin it is given and will stay spinning at 9 ~ 0. 

(ii) Group II (l--§ < £ < l + -§ ) 

Commercial tippe tops belong to Group II. A typical graph of vpy for a tippe top of 
Group II is shown in Fig. 12(b). We see that vpy is positive for 0<#<7r. Therefore, 
the end point at 9 = is unstable while the other end at 9 = tt is a stable point. 
Once given a sufficiently large spin at 9q 0, the tippe top of Group II will turn 
over and spin at 9 = it. Actually, commercial tippe tops have stems. Thus for those 
tops the above statement is valid up to the angle when the stem touches the table. 

(iii) Group III (1 + | < £ ) 

For an example of the tippe top of Group III, we may imagine a prolate spheroid 
put inside of a hollow massless sphere. The symmetric axis is along the diameter 
of the sphere and the mass distribution is nonuniform so that the center of mass is 
apart from the sphere's center. Figure |2](c) shows a typical graph of vpy for a tippe 
top of Group III. Similarly to the case of Group I the graph crosses the line vp = 
at an angle 

g / = C0S " 1 L (1 !_ A) l and \< Q f<^- (3- 11 ) 

In this case vpy is positive for < 9 < Of and negative for Of < < n. So, both 
ends at = and n are unstable points, while the angle Of is a fixed point. When the 
body is spun sufficiently rapidly with the initial angle anywhere, will approach 
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the fixed point 9f. Thus the tippe top of Group III, even though given a sufficiently 
large spin at 9q ~ 0, will never turn over to 9 = n but up to the angle 9f. 

Now it should be emphasized that the argument so far for the classification of 
tippe tops into three groups is based on the assumption that the GBC, £ = 0, is 
exactly satisfied. It is very interesting to note that the above classification into 
three groups and its classificatory criteria totally coincide with those obtained by 
Hugenholtz [3] and Leutwyler j^J, both of whom resorted to completely different 
arguments and methods. In fact, Hugenholtz considered the effect on the tippe 
top when a small frictional force is working during the uniform motion and reached 
the same conclusion. On the other hand, Leutwyler used Lagrangian formalism 
and searched for the minimum of energy for the tippe top under the constraint of 
Jellett's constant (|2.7|) . Finally the behavior of the tippe top under the GBC was 
studied earlier by Sakai Unfortunately, his work was written in Japanese and 
is, therefore, not well known. The consequences derived in this subsection partly 
overlap with his results. 

3.3 The behavior of the variable £ 

As stated before, the GBC is not satisfied initially for the tippe top. Actually, the 
initial value of £ is large and positive. We have performed numerical computations 
to see the behaviors of £ and 9 in time t. Typical examples are shown in Figs. Eland 
01 where the scale of the left sides is for £ normalized by the initial value £o> while 
the scale of the right sides is for 9 in radian. Input parameters are for both cases 

R = 1.5 cm, a = 0.15 cm, M — 15 g, g — 980 cm/sec 2 , 

A = C=-MR 2 , /i = 0.1, A = lcm/sec. (3.12) 
5 

For initial conditions we choose no = 100 rad/sec, 9q — — 0, and uq = for both 
cases, but we take 9 — 0.1 rad for the simulation shown in Fig|3] and 9 = 0.01 
rad for the one in FigEJ With these initial conditions we have £o = Cn$ and 
J = £o(-Rcos#o — a). We show in FigElthe trajectories of slip velocity of the contact 
point P in the (vpx, vpy) space which are obtained from the above simulations with 
(a) #o = 0.1 and (b) #o = 0.01. The argument in section 3.2 tells us that a tippe top 
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represented by the input parameters (|3.12|) is classified into Group II and, therefore, 
this tippe top would turn over up to an inverted position, 9 = n, when it is given a 
sufficiently large initial spin. 

Also plotted in Figs|3] and |U are the curves £ + and the expressions of which 
are given below. Using Eqs. (j3.2H3.3|) and replacing Q and L z with £, 9 and J, we find 
that the X- and Y- components of the rotational equations (j2.5aH2.5b"j) are rewritten 
as 

£sin# = U(£,9,J)0 + h(9)F Y , (3.13a) 
AO = - J ^V^,9,J)sm9-h(9)F x , (3.13b) 

with 

ttic a j\ J-£{Rcos9-a) 

U(£,0,J) = £cos#, (3.13c) 

V(€,9,J) = {J -£(R cos 9- a)} £ + AMgah{9) , (3.13d) 

where we have set N = Mg. We see in Figs. OlandEJ especially in the former, that 
9 changes while nutating. We also see that the inflection points of 9, where the 
condition 9 = is satisfied, fall on a rather smooth curve about which 9 nutates. 
At these inflection points of 9, the right-hand side (RHS) of Eq. (j3.13b|) vanishes. 
Here we note that, unless sin9 m 0, the second term —h(9)Fx may be neglected as 
compared with the first term, since Fig|S] shows the smallness of vpx- Then solving 
V(£, 9, J) = for £, we obtain 



J ± a/J 2 + 4AMg(R cos 9 — a) (R — a cos 9) a 
e± = 2(Rcos9-a) • (3 ' 14) 

We expect that at the inflection points of 9 and if not sin# ~ 0, £ takes the values 
which are either on the curve £ + or In the limit (AMgR 2 a) / J 2 <C 1 and a <C R, 
which is true in these simulations, we have 

u » 7r ^, e- » , (3.i5) 

Kcos9 — a J 

and thus £_/£o ~ 0. 

Fig 01 shows the result of the simulation with an initial value 9q = 0.1 rad. The 
variable £, starting from a large positive value £o = Cno, begins to fluctuate around 
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(i) (ii) 

Figure 3: (i) Time evolution of the variable £ and the inclination angle 9. Input 
parameters are R= 1.5 cm, a = 0.15 cm, M— 15 g, g = 980 cm/sec 2 , A = C = 2/5MR 2 , 
[1 = 0.1, A=l cm/sec. Initial conditions are 9 = 0.1 rad, n = 100 rad/sec, # = ^o = 0, 
Uq = 0. The curves £± are given by Eq. ()3.14|) . (ii) Blow-up of the section surrounded 
by dashed lines in (i). 
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Figure 4: (i) Time evolution of the variable £ and the inclination angle 6 with an 
initial condition 9q = 0.01 rad. Input parameters and other initial conditions are 
the same as in Figure El The curves £± are given by Eq. (j3.14|) . (ii) Blow-up of the 
section surrounded by dashed lines in (i). 
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Figure 5: The trajectories of slip velocity of the contact point P in the (vpx,vpy) 
space which are obtained from (a) the simulation shown in Fig|3]with 6 = 0.1 and 
(b) the one in FigHwith 9 = 0.01. 

the curve £+. The fluctuation of £ becomes larger but £ is still positive for a while. 
The inclination angle 9 decreases rapidly from the initial value 9 and then starts 
to nutate. The amplitude of nutation becomes larger and the minimum value of 9 
decreases further. And at a certain point where 9 ~ 0, the fluctuation of £ becomes 
so large that £ takes negative values. Then £ starts to fluctuate around the curve £_ 
and 9 is going to increase while nutating. The fluctuation of £ is getting smaller as 
9 is increasing, but it becomes large again when 9 approaches n. We have observed 
in Fig 01 that the fluctuation of £ around the curve £ + at the beginning soon shifts to 
the one around the curve For this rapid transition of £, the simulation shows that 
the system should pass through the phase where 9^0. When a simulation starts 
with a very small initial value 9$ as in Fig0J then £ quickly moves to a fluctuation 
around the curve 
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Let us look more closely the behavior of £ in Fig|3] at an early stage (to be 
specific, between 0<£<1 sec). Recall J = £o(i2cos0o — a). Then Eq. (j3.13cj) gives 
[/(£, 9, J)\t=o = — ^ocos^o, which is large and negative. At the very beginning of 
time, FigEl (a) shows that v PY /\v P \ « 1, and thus we have F Y ~ —\iMg. Also the 
term £ sin 9 on the left-hand side (LHS) of Eq. (J3.13aj) may be neglected in the leading 
order as compared with the h(9)Fy term, since (d(£,/£o)/dt)sm9 ~ sin6> x (1/sec) 
and sin 9 is small, while (iMgR/^Q ~2 x (1/sec). Hence we find from ()3.13a|) . 

9 ~ j Tt ^ n 1 < , at the very beginning, (3.16) 
U{€,0,J) 

which explains a rapid decrease of 9 from an initial value 9q. 

Along with the rapid decrease of 9, FigJ3] (a) shows that the slip velocity Vp 
of the contact point P tends to vanish. Then, in this region where 9 is small and 
v P « 0, the term h(9)F Y of the RHS of (J3.13aj) may be neglected while U(£,0, J) 
is expressed as U(£, 9, J) ~ (£o — 2£). Hence Eq. (|3.13ajl is reduced to 

£sin0 = (£ Q -2£)0, (3.17) 

and its solution is given by 

1 + cos 9 

2£ - £o = const, x . 3.18 

1 — cos 9 

We observe in the simulation shown in FigJSlthat the behavior of £ during the time 
0.05 < t < 0.4 sec is approximately described as 

f„l +Cl l±£^, (3,9) 
£ 2 1 — cos 9 

with a positive constant G\. Although the term h(9)Fy has been neglected to 
derive (J3.17)) . the small effect of the frictional force still remains and it produces 
the nutation of 9, which in turn gives £ a fluctuating behavior around the curve £+ 
according to ()3.19|) . Along with nutation, the minimum value of 9 further decreases 
and so the fluctuation of £ is getting larger. 

Then at a certain point (at t w 0.45 sec), the behavior of £ shifts to the one 
which is, later on up to 1 sec, roughly described as 

i »i-fti±^, (3.20, 
£o 2 1 — cos 9 
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with a positive constant C2, and £ may take negative values. Actually £ fluctuates 
rapidly between positive and negative values. Also, with the shift of the behavior 
of £, £/(£, 9, J) turns to always take positive values. In this region, vpy is small but 
positive on the average in time (see Fig El (a)). Now taking the time average of both 
sides of (|3.13aj) . we see 9, the time average of 9, is positive, since the LHS, £ sin^, 
may be neglected while h{9)Fy is negative. Thus, from 0.45 sec to 1 sec, 9 gradually 
increases while nutating. As 9 is increasing, the effect of sin# on £ in the LHS of 
()H.13a|) gets weaker and the fluctuation of £ becomes smaller. In the end £ oscillates 
mildly about a negative value 

When we start simulation with a very small initial value #0 = 0.01 rad as in 
FigEl £ quickly takes negative values and 9 starts to increase. The fluctuations of 
£ and 9 are much smaller than those in FigO With a smaller 9$, the center of 
mass O receives less recoil from the frictional force F, which explains the smaller 
fluctuations for £ and 9. 

So far we have shown the result of the numerical analysis for a tippe top which 
belongs to Group II. Given a sufficiently large spin, the GBC for this tippe top, 
which is not fulfilled initially, will soon be satisfied approximately, and the body 
will start to turn over. Actually the GBC, £ = 0, is modified to £ = £ m = —AMgaR/ J 
and |£ m /£o| 1. This modification has an only effect of shifting the positions of 9 C 
(EHUJ) and 9 f (pHT| slightly. 

Empirically we know that when a given spin is not fast enough, the tippe top 
does not turn over and stays spinning with its stem up. Later in Sec. 4, we argue 
that there exists a critical value for the initial spin given to the tippe tops of Group 
II and III. If the initial spin is below this critical value, then even the tippe tops of 
Group II and III do not turn over. We have performed similar simulations as those 
in Figs. El and 0] with the same tippe top and the same initial conditions, except 
that the initial spins are below the critical value. In these simulations we find that 
£, starting from a positive £ > fi rs t fluctuates around the negative value £_ and then 
returns to positive values and fluctuates around £+, while the inclination angle 9 
remains approximately zero. A typical example is shown in Fig. El where input 
parameters and initial conditions are the same as in Fig. El (and thus the tippe 
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Figure 6: Time evolution up to 12 sec of the variable £ and the inclination angle 9 
for a tippe top of Group II with an initial spin n =30 rad/sec. Input parameters 
and other initial conditions are the same as in Figure EJ The curves £± are given by 

Ea. icrm 

top for this simulation belongs to Group II), except that the initial spin n is 30 
rad/sec. The critical value for the initial spin is given by n\ in (j4.18j) below and we 
have rii =36 rad/sec for this case. 

We also performed simulations for the tippe tops of Group I, which are predicted 
to stay spinning at 9 « however large a spin they are given. Plotted in Fig. UJare 
the time evolution of £ and 9 for a tippe top belonging to Group I with initial spins 
(a) Uq = 100 rad/sec and (b) n = 30 rad/sec. Given a large initial spin (Fig. [7(a)), 
£ for a tippe top of Group I stays positive and takes values very close to £ + . But, 
with a small initial spin (Fig. UJ^b)), £ changes from a positive £o to negative and 
fluctuates around £_ for a while, and then back to positive values again. In both 
cases the tippe top stays spinning at 9 ~ 0. 

For completeness we show, in Fig. |HJ typical examples of the time evolution 
of £ and 9 for a tippe top belonging to Group III with initial spins, (a) n = 100 
rad/sec and (b) n = 15 rad/sec. Input parameters and other initial conditions are 
explained in the caption of Fig. |HJ The critical value for the initial spin for the tippe 
top of Group III is given again by n\ in ()4.18|) below and we have n\ =23.5 rad/sec 
for this simulation. When the initial spin no is larger than the critical value n\ (Fig. 
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Figure 7: Time evolution of the variable £ and the inclination angle 9 for a tippe top 
of Group I with initial spins (a) n = 100 rad/sec and (b) n = 30 rad/sec. Input 
parameters are R = 1.5 cm, a = 0.15 cm (and thus a/R = 0.1), M = 15 g, # = 980 
cm/sec 2 , A/C = 0.85, C = 2/5MR 2 , /i = 0.1, A = l cm/sec. Other initial conditions 
are 9 = 0.1 rad, 9 = Q = 0, u = 0. The curves £± are given by Eq. (j3.14| ). 



EJa)), the variable £ for a tippe top of Group III shows a similar behavior as the one 
presented in Fig. |3]for the tippe top of Group II. To be specific, £ becomes small and 
fluctuates around the curve £_ while 9 increases. Note that the tippe top of Group 
III never turns over to the inverted position, 9 = n. In the simulation of Fig. IHJa), 
9 goes up to the asymptotic angle flasympj which is below the fixed point 9 f (=2.21 
rad.) derived from (JXTTJ). (See also the discussion on the plot in Fig. ITUjI . Given a 
smaller initial spin than n\ (Fig. Efb)), £ for a tippe top of Group III, starting from 
a positive £o, fluctuates around the negative value £_ and then becomes positive 
and fluctuates around £ + , while 9 remains approximately zero, a similar behavior as 
the one shown in Fig. [HI for the case of a tippe top of Group II with an insufficient 
initial spin. 

From these numerical analyses we see that the behavior of £ is closely related to 
the inversion phenomenon of the tippe top. As the tippe top turns over, simulation 
shows that £ becomes small (in the sense |£/£o| ~ \£,(R~ a )l J\ ^ 1) an d takes values 
close to ~ £ m , which implies that the relation (|3.4|) is approximately satisfied. 
Conversely, when the relation (J3.4j) holds, it means that the center of mass of the 
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(b) 

Figure 8: Time evolution of the variable £ and the inclination angle 9 for a tippe 
top of Group III with initial spins (a) no = 100 rad/sec and (b) no = 15 rad/sec. 
Input parameters are R = 1.5 cm, a = 0.225 cm (and thus a/R = 0.15), M = 15 g, 
g = 980 cm/sec 2 , A/C = 1.25, C = 0.8 x (2/5)Mi? 2 , /i = 0.1, A = 1 cm/sec. Other 
initial conditions are # = 0.1 rad, #o — — 0, u = 0. The curves £± are given by 
Ea. rfTm 



tippe top goes up as Lz decreases. 



4 Stability and critical spin 

4.1 Steady states 

In Sec. 3. 2 we have studied the behaviors of the spinning tippe top when the gy- 
roscopic balance condition £ = is exactly satisfied. The situation corresponds to 
the one in which the tippe top is given an infinitely large initial spin. Actually the 
initial spin given to the tippe top is finite and we know empirically that a tippe top 
with a small spin is stable and does not turn over. We will now consider how large 
an initial spin should be for the tippe top to turn over. For that purpose we will 
study the steady states of the tippe top and examine their stability. 

Actually the steady states (or the asymptotic states) of the tippe top and their 
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stability were analyzed by Ebenfeld and Scheck [ES] [7j. They used the total energy 
of the spinning top as a Liapunov function. Then the steady states were found as so- 
lutions of constant energy. The stability or instability of these states was determined 
by examining whether the Liapunov function assumes a minimum or a maximum 
at these states under the constraint of Jellett's constant. The tippe top inversion 
was also analyzed recently by Bou-Rabee, Marsden and Romero [BMR] [9] as a 
dissipation-induced instability. BMR used the modified Maxwell-Bloch equations 
and an energy-momentum argument to determine the stability of the non-inverted 
and inverted states of the tippe top. 

Here we take a different approach to this problem. And we discuss the stability 
of the steady states in terms of the initial spin velocity n given at the non-inverted 
position 9 = 0. Recently, Moffatt, Shimomura and Branicki [MSB] made a linear 
stability analysis of the spinning motion of spheroids They identified the steady 
states, and then discussed their stability and found the critical angular velocity 
needed for the rise of the body. In order to find the steady states of the spinning 
tippe top, we adopt the method taken by MSB for the case of spheroids. But for 
the stability analysis of the steady states, we develop a new stability criterion which 
is different from the ones used by ES, BMR and MSB. 

Our approach to the stability problem of the tippe top is as follows. Once a 
steady state is known, the system is perturbed around the steady state. Particularly 
we focus our attention on the variable 9, which is perturbed to 



where 9 S is a value at the steady state and 59 is a small quantity. Using the equations 
of motion, we obtain, under the linear approximation, a first-order ODE for 59 of 
the following form: 



where n s and Q s are values taken at the steady state. Equation ()4.2)1 implies that the 
change of 59 is governed by the sign of the function H. If H is positive (negative), \59\ 
will increase (decrease) with time. Therefore, we conclude that when H is negative 
(positive), then the state is stable (unstable). This is the criterion for stability of 



e = e s + 59 , 



(4.1) 




(4.2) 
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the steady state, which we will use in this paper. The stability criterion in this 
work is derived from an intuitive analysis of the equations of motion. We check in 
Appendices C and D that they are consistent with those derived by ES and BMR 
which are based on mathematically rigorous methods. 

Superficially the above criterion (|4.2jl seems quite different from the one used 
by ES [7], but actually we have found that both are equivalent and, therefore, our 
results are consistent with theirs. In Appendix C we will show the equivalence of 
both criteria and that the stability conditions of the steady states which we will 
obtain coincide with the ones found by ES. After all, ES utilized the total energy 
(an integral form) [7j, while we will use equations of motion (differential forms). 

The criterion (J4.2j) for the stability of the tippe top also seems different from the 
ones used by BMR j^j, which were derived from the tippe top modified Maxwell- 
Bloch equations. In order to obtain the stability criteria, both BMR and we linearize 
equations of motion about the steady states and use sliding friction, which is assumed 
to be an analytic function of the slip velocity, as the main mechanism behind tippe 
top inversion. Thus it is well expected that both criteria lead to the consistent 
results on the stability of the non-inverted and inverted states. (The stability of the 
intermediate states have not been analyzed yet by means of the modified Maxwell- 
Bloch equations). In Appendix D we will show that the expressions of the criteria 
provided in BMR become more transparent when they are rewritten in terms of 
the parameters and classification criteria used in this paper, and that they lead to 
the same stability conditions for the vertical spinning states which will be obtained 
later by using the criterion (|4.2|) . Besides, although BMR did not mentioned, the 
classification of tippe tops into three groups, Group I, II, and III, is shown to be 
possible through the close examination of the criteria in BMR. 

The steady states of the spinning motion of the tippe top are obtained from the 
equations of motion (j2.6aM2.6c|) and (j2.11aH2.11c|) by setting Q = = 9 = n = iiox = 
uoy = uqz = [12]. Since we assume that the sliding friction (j2.14j) . i.e., a modified 
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version of Coulomb law 3 , is the only frictional force present, the energy equation 

^E = F . Vp = -ptN-^- (4.3) 
dt \vp(A)\ v 1 

shows that Vp = and F = at the steady states ^2] • Thus we obtain for the 
steady states of the tippe top, 

uox = 0, (4.4a) 

Qu OY = 0, (4.4b) 

Q(Cn - AQ cos 9) sin 9 + Mga sin 9 = 0, (4.4c) 

uoy + {R(n - n cos 9) + af2} sin 9 = 0, (4.4d) 

where = Mg and the velocity equations (j2.8aH2.8b|) and (|2.fJj) have been used. The 
solutions for Eqs. (j4.4aM.4dp are: 

i) Vertical spin state at 9 = : 

uox = Moy — , 9 = 0, n arbitrary, Q undefined , (4.5) 

which is a spinning state about the axis of symmetry with the center of mass below 
the sphere' center. 

ii) Vertical spin state at 9 = tt : 

uox = uoy = , 9 = n , n arbitrary, Q undefined , (4.6) 

which is an overturned spinning state about the axis of symmetry with the center 
of mass above the sphere' center. 



iii) Intermediate states: 



uox = uoy = , < 9 < n , 

tt(Cn - Ail cos 9) + aMg = , (4.7a) 
R{n - Q cos 9) + aVL = . (4.7b) 



3 Recall that the exact Coulomb friction (|2.13|) is non-analytic at vp = and a nonlinear friction 
law that would not appear in the linear approximation [S] . 
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The elimination of n from (|4.7a|) and (j4.7b|) gives 

= ^ . (4.8) 

(A-C)cos8 + C^ y } 

The necessary (but not sufficient) condition for the existence of such states is 

A \ „ CL , 

--1 cosfl + ->0. (4.9) 
G j ti 

Recall (jH.fjjl which was used for the classification of tippe tops into three groups in 
Sec. 3.2. Thus, intermediate states may exist at 8> c = cos -1 f ^^Sry l for the tippe 
top of Group I (~ < 1-f ), at 8 between and n for Group II < ^ < 

and at 8 < Qf = cos~ 1 ( ^fz^ j for Group III (1 + | < ±). 

There appear, in total, three categories of steady states for a loaded sphere 
version of the tippe top. 

4.2 Stability analysis of the steady states 

The turnover of the tippe top is associated with the effect of the sliding friction 
(with a coefficient /x) at the point of contact P. Near the steady states, we know 
that vp m 0, which is equivalent to the situation where \i ~ 0. Thus for the stability 
analysis of the steady states, we consider the limiting case of /i C 1 ^2]- Since we 
expect Jj~0(/i) near the steady states, we have 6~0(n) and 9~0(fj, 2 ). Eq. (|2.8b|) 
shows v TOt pY ~ 0(1), which leads to vpy ~ 0(1) and, hence, F Y ~ O(p). Then 
(J2.11bj) gives Uox~0(fJ'), and fpjsr~0(/i) from ()2.8a|) and (|2.9j) . and thus we have 
F x ~C(Ai 2 ) and u y~0(ij, 2 ). 

The above order estimation in /i near the steady states leads to the primary 
balance in (j2.6b|) which holds at leading order in \x [T2] . 

n(Cn - AVL cos 9) + aMg = . (4.10) 

Note that with a sufficiently large Q, Eq. (j4.10|) reduces to C, = Cn— AQ cos 9 = 0, the 
GBC. 



26 



4.2.1 Stability of the vertical spin state at 6 = 

The angle 6 is perturbed from 6 = 0, and we take 6 = 56 <C 1. In the linear 
approximation we may take n = const., since Eq. (|2.6c|) implies that h is quadratic 
in small quantities (note Fy ~ O(fi)). With cos 56 = 1, the primary balance ()4.10j) 
gives 

Q = — |Cn ± v 7 !^) 2 + 4AM#a} . (4.11) 
In this approximation f2 is also a constant. Then Eq. ([2.6aj) gives 

2AQ-Cn \v P (A)\ K ' 

and we may take 

VpY = {R( n - Q) + aVl} 56, (4.13) 
since uoy ~ 0{.^ 2 )- Hence, we require for the stability at 6 = 

Using the expressions of both "+" and "— " solutions for Q in (|4.11jl . the above 
condition is rewritten as 

± ( ^— — -C)n< ^(Cn) 2 + 4AM ga , (4.15) 
\R — a J 



which gives 

.{£_(!_£)}<**(!_£)". (4 , 6) 



n 2 



IC y R') C V 
It is easily seen that the requirement ()4.16|) is always satisfied for any spin velocity 
n by the tippe top of Group I (^ < 1 — -|). As for the tippe top of Group II or III 
with ^ > (l — the requirement ()4.16|) is rewritten as 



The stability of the vertical spin state at 6 = is summarized as follows: For 
the tippe top of Group I with ^ < (l — ~), the spinning state at 6 = is stable for 
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any spin n, while for the tippe top of Group II or III with ^ > (l — -|) we require 
n < rii for its stability. In other words, the tippe top of Group II or III becomes 
unstable at 9 = if 



Mga ( a 



n ( 9 = )>„ 1= y c{ ,_»_ f)} ll -- j (4.18) 
4.2.2 Stability of the vertical spin state at 9 = n 

A similar analysis can be made for the stability of the spinning state at 9 = n. Now 
put 9 = 7i — 89' with 69' <C 1. Again we may take n = const., but note that n may 
be negative near 9 = ix. With cos# = — 1, the primary balance (j4.1U|) gives 

Q, = — {-Cn ± ^{Cnf - AAMga} . (4.19) 
In order for Q to have a real solution, we require 



2JAMga , A , 

\n\ > n = n 2 . (4.20) 

For | n | < 7i2, the spin is insufficient to overcome the effect of gravity and the 
orientation becomes unstable [T2*] . 
With Q = 0, Eq. (gUj) gives 

M = - ^^^ ^Mgr^r , (4.21) 

2AQ + Cn \v P (A)\ v ' 

and we may take, 

VpY = {R( n + + aVl} 59'. (4.22) 



Hence we require for the stability at 9 = n, 

R(n + ft) + afi 



> . (4.23) 



2AQ + Cn 

Using the expressions of both "+" and "— " solutions for Q in ([4. 19)1 . the above 
condition gives 

n . {(1 + £ ) _^ }> ^( 1 + £)'. ( , 24) 
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First, the requirement (|4.24|) is never satisfied by the tippe top of Group III 
> 1 + So the tippe top of Group III is unstable at 9 = n. Actually it never 

turns over to the position with 9 = ir. For the tippe top of Group I or II which 

satisfies 4 < (1 + -|), the requirement (|4.24|) becomes 



c V- iJ> 

^> ^ t ,"Ts ^ (1 + 9) =n2- (4.25) 



2 Mga / a \ 2 

" > C{(l + ^)-^}l 1 + ^ 



Note that rig > n|. 



The stability of the vertical spin state at 9 = ir is summarized as follows: For 
the tippe top of Group III (1 + j| < ^), the spinning state at 9 = it is unstable for 
any spin n, while for the tippe top of Group I or II with ^ < (1 + -|), the state at 
9 = tt is stable if 



H9 = tt)| > J Mga A (l + ^)=n 3 . (4.26) 

4.2.3 Stability of the intermediate state 

We have learned in Sec. 4. 2.1 that the spinning state of Group I at 9 = is stable. 
We also know from the discussion in Sec. 4.1 that the intermediate steady states of 
Group I, if they exist, must occur at 9 > 9 C = cos -1 (^ R ^ a _A^ j- This implies that 
the spinning motion of Group I near 9 = does not shift to a possible intermediate 
steady state. On the other hand, the tippe tops of Group II and III become unstable 
at # = when they are spun with a sufficiently large initial spin n(9 = 0) > ni, where 
rii is given by ()4.18j) . and they will start to turn over. Here we are interested in 
the intermediate steady states of the tippe top which are reached from the initial 
spinning position near 9 = 0. Therefore, in this subsection, we focus on the possible 
steady states only for the tippe tops of Group II and III, and examine their stability. 
The Jellett's constant given by ()2.7j) or (jH.Hjl is rewritten as 

J = Cn(Rcos9-a) + AQR sin 2 9 . (4.27) 

Now Eqs. f|4.7a|) and ()4.7b|) and the above expression of J completely determine the 
intermediate steady states. They are derived by solving 



r (^-l)cos# + -^ = <|(cosfl-^) 2 + ^sin 2 #j> , (4.28 



C R 



a NO A 



R J C 
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where 



Define the following function: 



where x = cos 6 and 



J 2 

MgaCR 2 



(4.29) 



Fix) = ||) , (4.30) 



fi(x) = (^-^ + 1' ( 4 - 31a ) 
= {(x-|) 2 + ^(l-a; 2 )} 2 . (4.31b) 

Then, Eq. (J4.28|) is rewritten as 

F(x) = re . (4.32) 
Since /^(x) = — 4\/^(x)/i(a;), we obtain 

F'(x) = ^v^J-^^^^-l), (4.33a) 
^V) = |^^{([/iW] 2 + (^-l)v / ^)) 2 + 3[/ 1 (x)] 4 }>0. (4.33b) 

The condition for the initial spin n(6 = 0) > n\ means J > Cn\(R — a). Using 
(ETTHJ) . we find re > (l--|) 4 /(^-l + ^), which leads to re > So we are looking 

for solutions of F(x) = re with re > -F(l). 

(i) Group II (l-!<£<l + i) 

When 1 < ^ < (1 + -|) , -F^ 3 -) < an d -F(^) is a monotonically decreasing function 
for — 1 < x < 1 . Hence, there is one and only one solution of F(x) = re at x s 
between —1 and 1, provided F(l) < re < F(— 1) . Otherwise, there is no solution, 
which means that there exists no intermediate steady state. Expressing J with the 
initial spin at # = as J = Cn(9 = 0)(R — a), we find that the condition re < F(— 1) 
gives 



(l-i) 
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Thus, in the case 1 < ^ < (1 + -|), one intermediate steady state exists at x s , 
provided that 

m < n(0 = 0) < n 4 . (4.35) 

We know F"{x) > from (I4.33bjl . and so F(x) is concave upward for — 1 < x < 1. 
If (1 — J|) < ^ < 1 , then may have a local minimum at a certain x between 
— 1 and 1. Recall that we are looking for the steady states which are reached from 
the position near 8 = and that the requirement for this is k > F(l). Hence, for 
the existence of such a steady state we need 

F(-l) > F(l) and F(-l) > k > F(l) . (4.36) 

The first condition F(—l) > F(l) gives 

A >1 _« (l + f) 4 -(l-f) 4 =r (437) 

and the second one F( — l) > k > F(l) leads to n\ < n(8 = 0) < n 4 . Some tippe 
tops of Group II with ~ < 1 satisfy > as well as the conditions ()4.36IK 

and thus r c < ^ < 1. For such tippe tops, the corresponding F(x) has a local 
minimum between Xd and 1, where Xd is a solution of F(xd) = F(l). These tippe 
tops, therefore, have one intermediate steady state at x s between —1 and Xd when 
the condition n\ < n{6 = 0) < n 4 is satisfied. See the discussion of case (c) in FigEl 
For the tippe tops of Group II with (1 — -|) < ^ < r c , there exists no intermediate 
state. We will see later, in the discussion of case (d) in FigEl that these tippe tops 
will turn over to 9 = tt once given a spin n(8 = 0) > ni, since F(—l) < F(l) and, 
hence, n\ > n 4 for these tops. 

(ii) Group III + A < 4) 

Since fi(x) should be positive, the allowed region ofxisx/<x<l with Xf = 
• Eq. ()4.33a|) together with (^ — 1) >0 shows that F(x) is a monotonically 
decreasing function for Xf < x < 1. Note that F(x) positively diverges when x 
approaches Xf from larger x. Hence, once k > F(l) is satisfied, F(x) = k has one 
and only one solution at x s such that Xf < x s < 1 . In other words, one intermediate 
steady state always exists at 6 S (= cos -1 x s ) between and 0f(= cos _1 x/) for the 
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tippe top of Group III, if the condition n{9 — 0) > n\ is satisfied. When n(9 = 0) 
gets larger, the angle 9 S gets closer to Of but never crosses 9f. In order for 9 S to 
reach 9f, n(9 = 0) should be infinite. 

Now we know that there exists an intermediate steady state for the tippe top of 
Group II with property r c < ^ < when n{9 = 0) satisfies rii < n(9 = 0) < n 4 . Also 
there is an intermediate steady state for the tippe top of Group III with < ^ 

if n(9 = 0)>ni. Let (n s ,Q s ,9 s ) represent such a steady state so that (n s ,Q s ,9 s ) are 
related by Eqs. (|4.7a|) and (|4.7b|) . and suppose this state to be perturbed to 

n = n s + 5n , ft = ft s + <5ft , 9 = 6 S + 59 . (4.38) 

Noting that 9 S = and Fy\ s = 0, we find that the perturbed state satisfies 

xa ^ Mg R 2 sin 2 9 S nf Ufl 

50 = t p r-D(x s )<J0 , (4.39) 

A c[S* + (£sm9 s y} 

where 

D(x s ) = 4 [A^)] 2 + - l) V/ 2 (x s ) . (4.40) 

The details of the derivation of ()4.39|) are given in Appendix B. 

If ^ > 1, then D(x s ) > 0. Also when r c < ^ < 1, we find that D(x s ) is 
still positive (see Appendix B). Thus we observe from ()4.39j) that 59 oc 59 with a 
negative constant at the intermediate steady state, which means that this state is 
indeed stable. 

Finally it is emphasized that the spinning state of the tippe top of Group I is 
stable at 9 = and the top will not turn over from the position near 9 = 0. On the 
other hand, the tippe top of Group III, when given a sufficiently large spin near the 
position 9 = 0, will tend to turn over and approach the steady state at 9 S but never 
up to the inverted position at 9 = ir. 



4.3 Critical spin for inversion of the tippe top of Group II 

The tippe top of Group II will turn over to the inverted position at 9 = n when it 
is given a sufficient initial spin. Let us estimate the critical value n c so that the 
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spinning top with n(6 = 0) > n c reaches the inverted position. 4 Recall that Jellett's 
constant ()4.27|) is invariant during the turnover from 9 = to 9 = n. From the 
relation Cn(6 = 0)(R— a) = Cn(8 = ir)(—R—a), we obtain 

n(9 = tt) = -^f^n(e = 0) . (4.41) 
R + a 

We already know that we need \n(9 = n)\ > n 3 for the stability at 9 = n, where n 3 is 
given in ()4.26|) . Thus we find 



"(fi = °) > l/ ^ M f? a,\ "n = ■ ( 4 - 42 ) 



Also from the instability condition of the tippe top of Group II at 9 = 0, we need 
n(9 = 0) > ni, where ni is given by (j4.18j) . Hence the condition for the tippe top 
of Group II to turn over up to 9 = tt is that the initial spin n{9 = 0) should be 
larger than both and n\. In fact, we observe > n\ for the tippe top with 
r c < ^ < ! + -§) while < ri\ for the tippe top with 1 — -| < ^ < r c , where r c is 
given by ()4.37j) . Therefore, we obtain 

n c =h fOT rc< ^< 1 + ^ (4.43) 
\n u for l-|<^<r c . V ; 

4.4 Numerical analysis 

We now study the time evolution of the inclination angle 9 from a spinning position 

near 9 = 0. Simulations are made with various values of ^ and j|, changing the 

input parameters A and a. Other input parameters are the same as those given in 

(j3.12|) . Initial conditions are #o = 0.01 rad, #o = ^o = 0, and it = 0, and the initial 

value of the spin velocity Uq is varied. Since we have chosen a very small 9 , we may 

consider n Q as n(9 = 0). 

Figure |H1 shows the asymptotic (final) angle of inclination, 9^^, as a function 

of n for several types of tippe tops of Group II with different values of ^ and 

4 The idea is borrowed from Ref.J21> where MSB estimated the critical angular velocity above 
which a uniform prolate spheroid will rise to the vertical state under the assumption of the GBC 
and thus the existence of Jellett's constant. 
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30 40 50 60 70 

n„ [rad/sec] 

Figure 9: The asymptotic value #asymp as a function of the initial spin velocity n 
for tippe tops of Group II with various values of ^ and -|; (a) the one with ^ = 1 
and J| = 0.15; the others have ^ = 0.95 but different j| such as (b) -| = 0.15, (c) 
| = 0.125 and (d) | = 0.1. 

(a) the one with ^ = 1 and -| = 0.15; the others have ^ = 0.95 but different -| such 
as (b) -| = 0.15, (c) — = 0.125, and (d) — = 0.1. The asymptotic angle ^ asymp may 
be or ir, or 6 S , the angle of a possible intermediate steady state. 

The symbols •, o, o and x represent the results for the tippe tops (a), (b), (c) 
and (d), respectively, and the thin solid curves (a), (b) and (c) are the trajectories 
obtained by solving (J4.28|) . We observe that the numerical results fall on the pre- 
dicted curves. The values of ni(n 4 ), in units of rad/sec, for the tops (a), (b), (c) 
and (d) are 34.4(62.9), 42.1(54.5), 45.7(49.4) and 51.4(44.4), respectively. In each 
case we see that the spinning state near = is stable when n <n 1 . Once n gets 
larger than n 1; the state becomes unstable and the tippe top turns over up to the 
asymptotic angle # aS ymp- For the tippe tops (a) and (b) the values of ^ymp grow 
with no from to tt. On the other hand, the tippe top (c) satisfies r c < ^ < 1 with 
r c = 0.94, and thus the intermediate steady state exists only at 6> s (= 6> asymp ) with 
8 d <8 s < 7r, where 9 d is a solution of F(cos8 d ) = F(l). We find 6 d = 1.89. Thus 
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20 40 60 80 100 120 140 

n„ [rad/sec] 

Figure 10: The asymptotic value #asymp as a function of the initial spin velocity no 
for tippe tops of Group III: (a) with £ = 1.25 and | = 0.025; and (b) with ^ = 1.25 
and -| = 0.15. 

when n gets larger than n\ for the case of the tippe top (c), the asymptotic angle 
#asym P jumps from to When n > n 4 , ^asymp = tt for the tops (a), (b) and (c). 
In the case of the tippe top (d), we find r c = 0.96 and thus ^ < r c , which leads to 
ni>rii. Therefore, there is no intermediate steady state, and the asymptotic angle 
#asym P is or Ti depending on n «s n\. 

We plot in FigjTUlthe asymptotic angle ^asymp as a function of no for the tippe 
tops of Group III; (a) with £ = 1.25 and | = 0.025 and (b) with £ = 1.25 and 
= 0.15. The symbols • and o represent the results of simulation for the tippe tops 
(a) and (b), respectively, and the thin solid curves (a) and (b) are the trajectories 
obtained by solving ()4.28|) . We observe again that the numerical results on 6* asym p for 
both tops (a) and (b) fall on the predicted curves. The values of rii for the tops (a) 
and (b) are 13.3 and 23.5 rad/sec, respectively. In both cases the spinning position 
near 9 = is stable when no is below n\. Above ni, the value of # a symp grows with 
no and approaches the fixed point Of. The values of Of for the tops (a) and (b) are 
1.67 and 2.21 rad, respectively. 
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Figure 11: The time evolution of the angle 9 for a tippe top of Group II from a 
spinning position near 9 = 0. 

For simulations we have used a modified version of the Coulomb friction F given 
in (|2.14jl . The value ^ asymp is not affected by the strength of the coefficient \i. The 
strength of fi instead has an effect on the rate of rising of the tippe top. If we use 
another form than (|2.14|) for the sliding friction, and moreover, it is expressed as a 
continuous function of Vp and vanishes at Vp = 0, then we still expect that we get 
the same numerical results on ^asymp vs. uq as shown in FigJHland FigHTA This is 
due to the observation that the numerical value # a symp has fallen on the predicted 
curves which are derived from (|4.28J) and that we have obtained (I4.28J1 using the 
property of F which vanishes at the steady states together with vp. 

Figure ^2 shows the time evolution of the inclination angle 9 for a tippe top of 
Group II from a spinning position near 9 = for various values of the initial spin 
velocity no- Input parameters and initial conditions are the same as before and we 
take ^ = 1 and j| = 0.15. The asymptotic angles ^asymp which will be reached are 0, 
0.92, 1.67, 2.49, tc and 7r rad for n =30, 40, 50, 60, 70 and 80 rad/sec, respectively. 
Simulations with a modified version of the Coulomb friction (|2.14j) show that the 
larger value of no is given, the faster the rate of rising becomes. 
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5 Summary and Discussion 



We have examined an inversion phenomenon of the spinning tippe top, focusing 
our attention on its relevance to the gyroscopic balance condition (GBC), which 
was discovered by Moffatt and Shimomura in the study of the spinning motion of a 
hard-boiled egg. In order to analyze the GBC in detail for the case of the tippe top, 
we introduce a variable £ given by ()3.1|) so that £ = corresponds to the GBC, and 
study the behavior of £. Contrary to the case of the spinning egg, the GBC is not 
satisfied initially for the tippe top. The simulation shows that, starting from a large 
positive value £ 0) the variable £ for the tippe tops which rise, soon fluctuates around 
a negative but small value £ m such that |£ m /£o| ~ 0. Thus for these tippe tops, the 
GBC, though it is not fulfilled initially, will soon be satisfied approximately. Once £ 
fluctuates around the value £ m , these tops become unstable and start to turn over. 
On the other hand, in the case of the tippe tops which do not turn over, £ remains 
positive around £ or changes from positive £o to negative values and then back to 
positive values close to £o again. 

Under the GBC the governing equations for the tippe top are much simplified 
and, together with the geometry of the tippe top, we obtain a first-order ODE for 9 
in the following form [TUj (see (j3.5|) or ()3.7^ : 

f t =K9). (5.1) 

It is noted that this equation has a remarkable resemblance to the renormalization 
group (RG) equation for the effective coupling constant g, 

ft = m ' (5 ' 2) 
which appears in quantum field theories for critical phenomena [To! IT7] and high 
energy physics [THJ. Here in ([5.2)1 . t is expressed as t = InA with a dimensionless 
scale parameter A. Provided that (3(g) has a zero at g = g c , we find that, if /3'(g c ) < 0, 
then g(t) — > g c as t —>■ oo (A —>■ oo), and while if /3'(g c ) > 0, g(t) — > g c as t — » — oo 
(A — > 0). The limiting value g c of g(t) is known as the ultraviolet (infrared) fixed 
point in the former (latter) case. Similarity between the two equations, (|5.1|) and 
(15. 2|) . and the notion of the RG equation brought us to a consequence that tippe tops 
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are classified into three groups, depending on the values of ^ and -|. A resemblance 
of Eq. (j4.2j) to the RG equation also gave us a hint that Eq. ()4.2|) might serve as a 
criterion for stability of the steady state in Sec. 4. 

The criterion (4.2) is a first-order ODE for the (perturbed) inclination angle 
58, and the results derived from this criterion coincide with those by ES and BMR 
which are obtained by mathematically rigorous methods. The key ingredients in 
the process of arriving at this first-order ODE are the order estimation in fi near 
the steady states and an intuitive analysis of the equations of motion. The criterion 
(4.2) can also be applied to the stability analysis of other spinning objects. In fact 
we have applied (j4.2j) to the spinning motion of spheroids (prolate and oblate) which 
was recently examined in detail by MSB ^2], and we have obtained consistent results 
with theirs. 

Finally we have assumed, in the present work, a modified version of Coulomb 
law (j2.14|) for the sliding friction, since Coulomb friction (j2.13J) is non-analytic and 
undefined at vp = 0. On the other hand, Cohen used Coulomb friction in his 
pioneering work on the tippe top j3], and analyzed its spinning motion numerically 
for the first time. He reported the result of a sample simulation in Fig.5 of his 
paper The Coulomb friction is realistic provided that \vp\ is away from zero, 
but its application to the spinning motion of the tippe top is very delicate. Near 
steady states (i.e., near 9 = or 7r or 9 S ), Vp almost vanishes (see, for example, Fig|3] 
(a) and (b)). And there the X- and ^-components of Vp/\vp\ are changing signs 
rapidly and moreover non-analytically, and so are the components of friction, Fx 
and Fy. Coulomb friction may not be adequate to be applied to such a situation. 
In fact, Kane and Levinson argued against the work of Cohen, because it did 
not include adequate provisions for transitions from sliding to rolling and vice versa. 
They reanalyzed the simulation of Cohen, assuming Coulomb law for sliding friction, 
but also providing an algorithm that rolling begins when \vp\ <e (with e <C lm/sec) 
is satisfied, together with another algorithm for the transition from rolling to sliding. 
They found that a transition from sliding to rolling occurs soon after the motion 
has begun and that values of 9 remain below 0.077 rad thereafter. Or [5] adopted 
a hybrid friction law adding viscous friction, which is linearly related to Vp, to 
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Coulomb friction. Other frictional forces such as the one which is due to pure 
rotation about the normal at the point of contact might have some effect. After all 
it is safe to say that we have understood general features of the tippe top inversion. 
But it would be not until we have had thorough knowledge of frictional force that 
we completely understood the inversion phenomena of the tippe top. And yet, it 
flips over. 
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Appendix 

A Equations of motion for the tippe top 

We enumerate the equations of motion which are used to analyze the spinning motion 



of the tippe top: 

A(l sine = (Cn -2 AQ cos 9)9 + (a - R cos 9) F Y , (A.l) 

AO = -n(Cn- Aft cos 6) sin6 -a sin 6N -h(6)F x , (A.2) 

Ch = Rsin6F Y . (A. 3) 

Miiox = MVtuov + F x , (A.4) 

Miiov = -MVtuox + F Y , (A.5) 

Miioz = N -Mg . (A.6) 



F = -^Ny^f—, with \v P (A)\ = Jv 2 PX + v PY + A2 ( A . 7 ) 

\Vp{A)\ v 

v PX = u ox -h(9)9, (A.8) 

vpy = uoy + {R(n — VL cos 9) + aVl} sin 9 , (A. 9) 

h(9) = R-acos9 (A.10) 

u z = asin6 9 (A. 11) 



B Stability of the intermediate state 

There exists an intermediate steady state for the tippe top of Group II with property 
r c < ^ < l + when an initial spin n(9 — 0) satisfies n\ < n(9 — 0) < n A . There is 
also an intermediate steady state for the tippe top of Group III if n(9 — 0) > n x . In 
this appendix we show that these steady states are stable. 
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Near the steady states the primary balance condition (|4.1(J|) holds at leading 
order in \i. Differentiating both sides of (j4.10j) with respect to t, we obtain 

(Cn - 2AVL cos 6)(l + CVLh + Att 2 sin 66 = 0. (B.l) 

Using (|2.6a|) and (|2.6c|) . and eliminating fl and n, we find 

-F Y \ (a- Rcos9)(Cn-2AQ cos 6) + AQRsm 2 d\ 

(Cn - 2 Att cos 9) 2 + (AQsm9) 2 ' 1 ' ' 

Let (n s ,Q s ,9 s ) represent an intermediate steady state so that n s , Q s and 9 S are 
related by (|4.7aj) and (|4.7b|l . and suppose this state to be perturbed to 

n = n s + 5n , = tt s + 5tt , 9 = 9 S + 59 . (B.3) 

Since 6 S = and Fy\ s = 0, the perturbed state satisfies 

59 = -5F Y ™± f , (B.4) 



where 



Cn s -2Att s cos9 s A a 
S s = g~ = 2— cos# s - (cos 6 S - — ) , (B.5) 

T s = S s (cos6 s - ^) + ^ sin 2 6 S . (B.6) 

At leading order in /1, we have vpy = v TOt pY (recall Moy ~(9( / u 2 )), and thus we obtain 
from (l2~8b"J) . 

5r Y = ——5v TOtPY 



uMg 
A 



Rsm9 s {8n + (-| - cos 9 s )8n + fi s sin 9 8 5o} (B.7) 



Now we expect that the perturbed state still satisfies the primary balance con- 
dition (j4.1U|) . since A9 and Fx are 0(fi 2 ). Then a variation around the steady state 
gives 

5n - S S 5Q + -fl s sin 9 S 59 = . (B.8) 
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where S s is given by (jB.5|) . Also taking a variation of Jellett' constant (|4.27|) around 
the steady state (and then, of course, we have 5 J = 0), we obtain 



a A 

[cos 9 S — —)Sn + — sin 2 9 S 5Q + S S Q S sin 9 S 59 = 
R C 



From (jB.8|) and (|B.9|) . Sn and 5Q, are expressed in terms of 59 as 



8n 
5tt 



1 



T 

-L S 
1 



.4 



— <Sf + (-sm6 



C 



Q s sin 9 S S9 , 



A a 
~\ $* - — (costf, - —) sin 9,59 . 



Inserting these expressions into ()B.7|) . and then we obtain from (jB.4|) 

^Mg R 2 sin 2 9 S 



(B.9) 

(B.lOa) 
(B.lOb) 



59 



A 



C<| ) S s 2 + (^sin^)2 



D(x s )56 



(B.11) 



where 



D(x a ) =4 ,A(x s 



+ (^-i)v / iS) 

and x s = cos# s , and Eqs. (j4.31a|) and (j4.31b|) have been used. 

If 1 < ~, then D(x s ) > 0. Also when r c < ^ < 1, D(x s ) is still positive, which 
is explained as follows: The expression of ()4.33a|) shows that the function D(x) is 
related to F'(x) as 



(B.12) 



F'(x) 



-D(x) 



(B.13) 



When r c < ^ < 1, an intermediate steady state at x = x s exists provided F(— 1) > 
F(l) and F(— 1) > k > F(l) . At that point F'(x s ) is negative, and thus D(x s ) is 
positive. 
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C Equivalence between the criterion of ES |7j and 
Eq.fllH) 

Ebenfeld and Scheck [7] analyzed the stability of the spinning tippe top using the 
total energy as a Liapunov function and gave the stability criteria for the steady 
states. We take a different approach to this stability problem. First the system is 
perturbed around the steady state. Then, using the equations of motion and under 
the linear approximation, we obtain a first-order ODE for 59 of the form given in 
(|4.2jl . We make use of this equation and give a different stability criterion. In this 
appendix we show that both approaches are equivalent and thus they lead to the 
same conclusions on the stability conditions of the steady states. 

ES wrote the total energy of the spinning top as the sum of two terms (ES-(33)) 5 

E = EU(r )3 ,L l] ) + EW(f},L ± ,s li2 ) , (C.l) 

the second of which contains all the terms that will vanish at the steady states, while 
the first depends on rj^ = cos# and Jellett's constant J. In terms of the parameters 
used in this paper, and E^ are expressed as follows: 

Eil) = 2A^M +M ^~>- < C2 > 
B< 2 » = 1 -M(ul x +ul, Y + ul z ) + 1 -A^ 

c 



2C(1 -» 2 r RG(Vs) J ' 1 ' 



with 

G(Vs) = l~vl + j(Vb ~ |) 2 , Vs = cos9 (C.4) 

Note that ES set R=l. The condition dE^(9)/d9 = together with uox =uoy = 

uoz = leads to the three solutions of the steady states: (i) vertical spin state at 9 = 

()4.5|) . (ii) vertical spin state at 9 = n ()4.6|) . and (iii) intermediate states ()4.28|) . or 

equivalently, (j4.7aH4.7b"|) . It is recalled that we have obtained these solutions starting 

5 From now on, we write the equation of Ref. 7 as ES-(**). The Jellett constant A defined 
by ES is related to our J as A = J/R. 
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from equations of motion. At these steady states vanishes. For intermediate 
steady states, the factor {£ + J[i] 3 — ^(^3 — j^)]/ RG(r] 3 )} in (jC.3|) reduces to zero, 
due to (|4.7aH4.7E|) and Jellett's constant given in f)4.27j) . 

Now we show that the criterion, Eq. (|4.2)l . for the stability of the steady states 
is equivalent to the one derived by ES [Zj. For the stability analysis of the steady 
states, the order estimation in /i near the steady states is important, which has been 
pointed out by MSB in their work on the linear stability analysis of the spinning 
motion of spheroids [T2*] . As explained at the beginning of Sec. 4.2, near the steady 
states we have 4 ~ O(fi), vpx ~ CKz-O, vpy ~ v rot pY ~ C(l)> an d uoy ~ CK^ 2 )- 
Since E^ is already O(fx) (recall that it vanishes at the steady states), we have 
~C(/i 2 ), while Thus near the steady states, the energy equation 

()4.3j) is written at leading order in /i as 

Suppose the steady states to be perturbed to n = n s + 5n, Q = Q S + 5Q, 9 = 9 s + 59. 
Since 8 S = 0, we have 8 = 58, and is expanded as 

dE^ d 2 E^ 



5d + o({58) 2 ^, (C.6) 



d6 d8 2 

where we have used the fact ^jf-\ e=e =0. 

Meanwhile v ro tpy is shown to be expressed as 

VrotPY = v(n s ,n s ,e s )se . (c.7) 

Actually we have already obtained the expressions ()4.13|) and ()4.22|) for f ro tpy(~ 
vpy) near the steady states at # = and 8 = 71, respectively. Also near the interme- 
diate steady states, 5n and 5fl are expressed in terms of 58 as (JB.lOa)) and (JB.lObj) . 
respectively, and thus we obtain (jC.7|) . Now using Eqs. (jC.5|) - flC.7|) we find 

™ w V{n S) VL S) 8 s ) 2 1 

59 = ~^ Mg \vpjh)\ d2E(1) \ ' ( ) 

dO 1 0=0 s 

which means that we can identify H in ()4.2j) as 



TT ( „ . A a \ a negative constant . 
H[n s ,n s ,8 s ,-,-j = . (C.9) 

' de 2 \e=e a 
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Hence we conclude that the following assertions are equivalent: a steady state is 
stable (unstable) <^=> H is negative (positive) <^=> ^p2~ \ e _ e is positive (negative). 

In fact, ES showed that if the quantity (ES-(39)) with the upper sign is positive, 
then |g_g is positive at 6 S = and the non-inverted rotating motion is Liapunov 
stable. On the other hand, starting from the equations of motion we derived H and 
obtained the condition (j4.16|) for the stability of the rotating motion at 9 S = 0. It 
is easily seen that the statement that the quantity (ES-(39)) with the upper sign is 
positive is equivalent to the inequality given in (|4.1fjjl . once we know that Jellett's 
constant at 9 S = is given by J = Cn(Q s = 0)(R — a). Similarly, if the quantity 
(ES-(39)) with the lower sign is positive, then ^py- | e _ e is positive at 6 s = ir and 
the completely inverted rotating motion is Liapunov stable. The condition that the 
quantity (ES-(39)) with the lower sign is positive is equivalent to the inequality 
given in (j4.24|) . Note, this time, J = Cn(6 s = Tx)(-R-a). 

As for the intermediate steady state (—1 < cos 9 S < 1), ES stated that if the steady 
state exists and the quantity (ES-(40)) is negative, then \ e _ e is positive and 
the state is Liapunov stable. In Sec. 4. 2. 3 we have shown that the stability of the 
intermediate steady state is determined by the sign of D(x s ) given in ()4.40|) . Now 
it is interesting to note that D(x s ) is related to (ES-(40)) as follows: 



A 2 + 3[(A-C)x s + C^ 



2 



D(x s ) = ^ — — — x (ES.(40)) . (CIO) 



AC 2 

Hence the condition that the quantity (ES-(40)) is negative is equivalent to D(x s ) > 
0. We have seen in Sec. 4. 2. 3 that there exists an intermediate steady state for the 
tippe top of Group II and also of Group III. (We have not considered a possible 
intermediate steady state for Group I, since such a state, even if it exists, cannot 
be reached from the initial spinning position near 6 = 0.) For these steady states, 
we have shown, in Appendix B, that D(x s ) is positive and, therefore, the states are 
stable. 
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D Modified Maxwell- B loch equations and stabil- 
ity criteria 



Recently Bou-Rabee, Marsden and Romero [BMR] treated tippe top inversion as 
a dissipation-induced instability. They showed that the modified Maxwell-Bloch 
(mMB) equations are a normal form for tippe top inversion and, using the mMB 
equations and an energy-momentum argument, they gave criteria for the stability 
on the non- inverted and inverted states of the tippe top 0. Although we have 
not explored the connections between the mMB equations and the first-order ODE 
(|4.2jl for 89, we show in Appendix D that our results on the stability of the vertical 
spin states are consistent with the criteria provided by BMR. Actually, rewritten in 
terms of dimensional parameters and classification criteria used in this paper, the 
expressions of those criteria become more transparent and they lead to the same 
stability conditions as ours for the vertical spinning states. Besides, although BMR 
did not mentioned, the classification of tippe tops into three groups, Group I, II, 
and III, according to the behaviors of spinning motion, is possible from the close 
examination of those criteria. 

BMR used the moments of inertia defined as the ones about the principal axes 
attached to the center of sphere instead of the center of mass. The correspondence 
between the parameters used by BMR and ones in this paper are as follows: 

a 1 - /ie* 2 A 



11 R ' a C 
7q"bmr — ' 11 



(D.l) 



J //|e*|Fr _1 2 Mga 

RC ' a BMR "T" 



where ^bmr is the spin rate of the initially standing equilibrium solution (we added 
a subscript BMR to distinguish from our Q), and the dimensionless BMR's "Jellett" 
constant, jq, is restricted to have a certain value, i.e., 7q = —(1 + e*). Also BMR 
expressed the vector from the center of sphere to the center of mass 5*0 (in the 
BMR notation OC) as SO = Re*k, where k is a unit vector along the symmetry 
axis. Using the tippe top modified Maxwell-Bloch equations, BMR obtained the 
stability criteria for the non-inverted state which are given by the three inequalities 
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in (BMR-(5.3)) 6 . 

They took k = ez (upward) in (BMR-(5.3)). Since the non-inverted state has 
the center of mass below the center of sphere, we have e* = — |e*| = — -|, and thus 
Obmr = n(9 = 0). The first inequality of (BMR-(5.3)) is rewritten as ^ > 0, which 
is always satisfied. Apart from some irrelevant positive constants, the second and 
third inequalities are expressed, respectively, as 

Mga ( a\A / a\ 5 v 2 A 



[n(6 = 0)] 2 C 



/ a \ A / a\°u A / a\ .„ , 

( 1 -l)o + ( 1 -5)^-c + ( 1 -l) >0 - (D ' 2) 

A / a\ Mqa / a\ 2 i . . 

From these inequalities, we find: 

(ai) In the case ^ < (1 — -|), i.e., for the tippe top of Group I, the above inequalities 
are always satisfied. In other words, the non-inverted states (8 = 0) of Group 
I are always stable. 

(aii) In the case ^ > (1 — i.e., for the tippe tops of Group II or III, the inequality 



f)D.3j) is satisfied if 

W> = of< T .^ P --) • <D-n 



|2 ^ Mga / _ ay 

C{$ RJ 



which is the same requirement given in (j4.17|) for the stability of the tippe top 
of Group II or III. Note that the inequality (|D.2|) is automatically satisfied 
when both ^ > (1 — ~) and inequality (|D.3|) hold. 

Thus, the BMR criteria (BMR-(5.3)) lead to the same result as ours on the stability 
of the vertical spin state at 9 = 0. 

The inequalities (BMR-(5.3)), which were derived as the stability criteria for the 
non-inverted state, can also be used for the stability criteria for the inverted state, 
but with some replacements. Since k = ez (upward), the inverted state has the center 
of mass above the center of sphere. Thus we have e* = -| and Qbmr = —n(9 = 7r). 
6 From now on, we write the equation (**) of Ref. 9 as BMR-(**). 
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Changing variables in inequalities (jD.2|) and (|D.3J) as a — > —a, -| — > — -|, and 

[n(6 = 0)] 2 -> [n(^ = 7r)] 2 , we obtain 

!+ob+ 1 +oh-7^+ O+o >°- ( D - 5 ) 



n(0 = tt)] 2 <7 V RJC \ RJ a 2 C \ R 

for the stability for the inverted state. From the above two inequalities, we see: 

(bi) In the case ^ > (1 + -|), i.e., for the tippe top of Group III, the inequality 
(jD.6|) is never satisfied. Therefore, the inverted states (6 = ir) of Group III are 
always unstable. 

(bii) In the case ^ < (1 + -|), i.e., for the tippe top of Group I or II, the inequality 
(|D.6j) is satisfied if 

Mga ( a\ 



W = '>''> g «i + i)-3 } l 1 + 5J- (D7 > 

which is the same requirement given in 1)4.25)1 for the stability of the tippe top 
of Group I or II at 9 = n. The inequality (JD.5)) is automatically satisfied when 



both £ < (1 + |) and inequality (ITX61) hold. 

Thus, the BMR criteria (BMR-(5.3)) also lead to the same result as ours on the 
stability of the vertical spin state at 8 = it. 

Actually, BMR derived also the stability criteria for the inverted state, taking 
k = —ez, which are given by the three inequalities in (BMR- (5. 4) ) 7 . Of course, 
we can use them to obtain the stability conditions for the inverted state. Taking 
now Obmr = — T^ n (@ = 7r ) an d e* = — j| in the second and third inequality in 
(BMR-(5.4)), we reach the same conclusions, (bi) and (bii). 

BMR discussed in Ref. |0] about the heteroclinic connection between the non- 
inverted and inverted states of the tippe top. They used an energy-momentum 
argument to determine the asymptotic states of the tippe top and obtained the 



7 The second inequality should read as er(l + e*) 2 [cr(l - e*) - (1 - ^e* 2 )} + v 2 {\ - e*) 7 + (1 - 
e*) 3 ne*Fr~ 1 (1 — fie* 2 ) > 0. The error is traced back to the missing factor of (7°^°) in the 
expression of F in BMR-(4.2). 



48 



explicit criteria for the existence of a heteroclinic connection, which are given in 
Theorem 6.2 and the appendix of Ref. PJ. In terms of the classification criteria 
and conditions obtained in this paper, the statement in BMR on the existence of a 
heteroclinic connection can be restated as follows: (i) A tippe top must belong to 
Group II in order to have a heteroclinic connection, (ii) Further more, the initial 
spin n(9 = 0) should be larger than n\ (Eq. (|4.18Jl ) and n 4 (Eq. (|4.42Jl ) so that a tippe 
top becomes unstable at # = and reaches the inverted position. The requirements 
n(9 = 0) > n\ and n(9 = 0) > 714, respectively, correspond to the criteria r > and 
r 4 >0 in Theorem 6.2 in BMR. 
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